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Introduction 

Since  the  detonation  performance  of  an  energetic  material  is  directly  related  to  its  density 
[1],  the  latter  property  is  an  important  criterion  in  evaluating  prospective  energetic  compounds.  To 
the  best  of  our  knowledge,  the  crystalline  forms  of  1  -  3  have  the  highest  densities  that  have  been 
measured  for  any  unsubstituted  C,  N,  H  compounds:  1:  1.69  g/cm^  [2];  2:  1.715  g/cm^  [3];  3: 
1.674  g/cm^  [4].  It  is  accordingly  interesting  to  consider  these  as  possible  starting  pomts  for 
ener<^etic  molecules  such  as  4,  the  dinitro  derivative  of  3.  In  this  context,  a  potential  concern  in  the 
cases  of  2  and  3  would  be  any  tendency  to  rearrange  to  the  azides  5  and  6  [5-12].  This  point  shall 
be  addressed.  However  1  -  3  do  exist  and  form  stable  crystalline  materials.  It  should  be  useful  to 
investigate  some  of  the  properties  of  these  molecules,  and  simultaneously  to  try  to  gam  some 
insight  into  the  reasons  for  their  high  crystal  densities.  With  these  objectives,  we  have  computed 
and  analyzed  their  molecular  electrostatic  potentials,  as  well  as  those  of  4  -  6. 


4 


5 


6 


Methods  and  Procedure 

The  electrostatic  potential  V(r)  created  at  any  point  r  by  the  nuclei  and  electrons  of  a 
molecule  is  defined  rigorously  by  eq.  (1). 


V(r)  =  X 

A 


fpOO^ 

J 


(1) 


Za.  is  the  charge  on  nucleus  A,  located  at  Ra.,  p(r)  is  the  electronic  density,  which  we  obtain 
from  the  molecular  wave  function.  V(r)  is  a  real  physical  property,  which  can  be  determined 
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experimentally,  by  diffraction  methods,  as  well  as  computationally  [13],  It  has  proven  to  be  a 
useful  tool  for  predicting  and  interpreting  molecular  interactive  behavior  [13-15],  Regions  of 
positive  and  negative  potential  on  the  molecular  surface,  defined  as  the  0.001  electron/bohr^ 
contour  of  p(r)  [16],  are  indicative  of  sites  attractive  to  nucleophiles  and  electrophiles, 

respectively. 

In  recent  years,  we  have  developed  more  quantitative  methods  for  using  the  electrostatic 
potential  to  analyze  non-covalent  molecular  interactions  [17-21].  We  have  introduced  several 
statistically-based  indices  that  are  defined  in  terms  of  the  surface  potential,  two  of  which  are 
presented  in  eqs.  (2)  and  (3). 


n  =  -i|V(r;)-V, 


(2) 


(3) 


In  these  equations,  V(ri)  is  the  value  of  V(r)  at  point  ri  on  the  surface,  and  Vg  is  the  average  value 
of  the  potential  on  the  surface.  In  a  similar  fashion,  V+Ir,)  and  V“(r,)  are  the  positive  and  negative 

-a.  --  ^ 

values  of  V(r)  on  the  surface,  and  V  *  and  V ^  are  their  averages:  V ^  -  I  V+(rj)  and 

V-  =  ]iv-(rj). 

^  "  j=l 

n,  given  by  eq.  (2),  is  the  average  deviation  of  the  electrostatic  potential  on  the  molecular 
surface.  We  have  shown  that  11  is  an  effective  measure  of  local  polarity  (or  charge  separation), 
which  may  be  quite  significant  even  in  a  molecule  having  zero  dipole  moment,  e.g.  CO2  andp- 
dinitrobenzene  [17].  The  total  variance,  given  by  eq.  (3),  is  a  measure  of  the  spread  of  the 

values  of  the  surface  potential,  and  is  particularly  sensitive  to  positive  and  negative  extrema 
[18,21], 

In  this  work,  molecular  geometries  and  energies  were  computed  with  a  nonlocal  density 
functional  procedure  (GAUSSIAN  92/DFT  [22],  Becke  exchange  and  Lee,  Yang  and  Parr 
correlation  functionals  [23,24],  6-3 IG**  basis  set).  Using  these  structures,  the  electrostatic 
potentials  and  related  quantities  were  calculated  at  the  ab  initio  HF/STO-5G  level,  in  order  that  they 
be  compatible  with  such  data  obtained  earlier  for  other  molecules  [17,20,18,19,21,25].  We  have 
shown  that  the  resulting  trends  in  the  surface  electrostatic  potentials  are  the  same  when  the  HF/6- 
31G*  basis  set  is  used  [26]. 
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Results 

Energies  and  structural  data  are  shown  in  Table  1.  There  is  overall  satisfactory  agreement 
between  the  calculated  (gas  phase)  and  experimentally-determined  (crystalline)  bond  lengths  for  1  - 
3.  An  interesting  feature  of  1  is  the  near-constancy  of  the  bond  lengths  around  its  perimeter, 

contrasted  to  the  considerably  longer  internal  C-N  bonds. 

The  surface  electrostatic  potentials  of  1  -  6  are  shown  in  Figures  1-5.  Their  dominant 
features  are  the  extensive  positive  regions  above  the  heterocyclic  rings.  (Note,  in  Table  2,  the 
values  of  Vs,max>  the  most  positive  surface  potentials.)  Formally,  all  of  these  molecules  are 
unsaturated  and  thus  are  expected  to  have  pi  electrons  above  the  rings;  in  the  case  of  benzene,  for 
example,  these  produce  a  negative  potential  of  -10.5  kcal/mole  above  the  ring  [27].  However 
analogous  negative  regions  are  not  found  in  1  -  6,  because  of  the  effect  of  the  strongly  electron- 
attracting  ring  nitrogens;  this  has  already  been  observed  earlier  for  azines  (e.g.  7  - 10)  [27,28]. 
The  electron-withdrawing  NO2  groups  in  4  cause  the  ring  system  to  become  even  more  positive 
(Figure  3b  and  Table  2). 


8  9  10 


The  ring  nitrogens  on  the  peripheries  of  1  -  6  have  negative  potentials  of  varying  strengths 
associated  with  their  lone  pairs.  In  Tables  2  and  3  are  listed  their  most  negative  values  on  the 
molecular  surfaces,  Vs,min(N),  and  also  in  three-dimensions,  Vmin  (N).  The  triply-coordinated 
nitrogens  in  the  central  portions  of  the  molecules  have  no  negadve  potentials  at  all.  The  peripheral 
nitrogens  accordingly  have  some  basic  character  and  some  hydrogen-bond-accepting  tendencies, 
although  relatively  weak;  among  the  known  compounds  1  -  3,  the  strongest  in  both  respects 
should  be  1,  for  which  the  present  results  combined  with  our  earlier  correlations  [29,30]  suggest  a 
pKa  of  about  0.0  and  a  hydrogen-bonding  ability  similar  to  that  of  pyrazine,  7.  It  might  have  been 
anticipated  that  1  will  be  similar  to  r-triazine  (8)  in  these  properties,  but  the  negative  potentials  in  1 
are  actually  stronger  than  in  8;  e.g.  VminlN)  is  -80.0  kcal/mole  in  the  former  and  -74.6  kcal/mole 
in  the  latter  [28].  This  is  presumably  because  the  peripheral  nitrogens  in  1  withdraw  some 
electronic  charge  from  the  one  in  the  center  of  the  molecule,  which  is  formally  positive.  The  least 
basic  ring  nitrogens,  overall,  are  those  in  4,  which  have  lost  a  portion  of  their  electronic  charge  to 

the  nitro  groups. 

Table  2  also  contains  the  calculated  areas  and  volumes  of  1  -  4,  based  on  the  0.001  au 
contour  of  the  electronic  density,  as  well  as  several  properties  computed  from  the  surface 
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electrostatic  potentials;  U,  a",  a“  and  The  n  values  of  1  -  4  show  that  these  molecules  have 
high  degrees  of  internal  charge  separation,  as  can  also  be  seen  in  Figures  1-3.  For  1  -  4,  fl 
ranges  from  21  to  25  kcal/mole,  although  we  have  usually  found  it  to  be  between  2  and  15 
kcal/mole  for  organic  molecules  [17].  1-4  are  good  examples  of  how  local  polarity,  as  measured 
by  n,  may  differ  from  net  overall  polarity,  given  by  the  dipole  moment.  1  and  2  have  high  values 
of  n,  even  though  their  dipole  moments  are  zero,  due  to  symmetry.  Both  3  and  4  have  non-zero 
dipole  moments,  but  that  of  3  is  much  larger,  6.12  vs.  1.24  Debyes  at  the  HF/STO-5G//HF/3-21G 
level;  this  is  because  the  polarity  of  the  NO2  groups  in  4  is  opposite  in  direction  to  that  of  the  ring 
nitrogens.  However  the  H  values  of  3  and  4  are  nearly  the  same  (Table  2). 

The  magnitudes  of  are  also  large,  between  238  and  354  (kcal/mole)2  for  1  -  3,  and  423 

(kcal/mole)2  for  4,  whereas  we  usually  find  it  to  be  below  180  (kcal/mole)^  [18,19].  The  high 
values  of  for  1  -  4  suggest  strong  electrostatic  interaction  tendencies.  Although  Figures  1  -  3a 

9 

show  that  the  molecular  surfaces  of  1  -  3  are  largely  positive,  the  approximate  similarity  of  a;  and 
in  each  case  indicates  that  the  variation  in  the  positive  and  negative  potentials  is  roughly  the 
same.  In  marked  contrast,  the  positive  potentials  in  4  are  clearly  much  greater  in  magnitude  than 

9  9 

are  the  negative,  as  well  as  more  extensive;  this  is  seen  in  Figure  3b  and  in  the  fact  that  a"  »  a". 


Discussion 

The  possibihty  of  tetrazolo-azido  tautomerism  is  well-established  for  nitrogen  heterocyclic 
systems  [5-12];  for  2  and  3  it  is  represented  by  eqs.  (4)  and  (5): 


The  general  experience  has  been  that  the  tetrazolo  tautomers  appear  to  be  favored  by  a  high 
electronic  density  associated  with  the  central  heterocychc  ring  and  in  particular  with  the  nitrogen 
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that  is  common  to  both  rings  in  the  tetrazoio  form  [5-9];  thus,  electron-withdrawing  substituents 
tend  to  promote  the  azido  tautomer.  However  environmental  factors  can  be  extremely  important  [5- 
12];  for  example,  the  equilibrium  shown  in  eq.  (6)  has  been  observed  in  chloroform  and  in 
trifluoroacetic  acid  solutions  [6],  but  only  11  is  present  in  the  solid  phase  [31]. 


(6) 


From  our  computed  energies  in  Table  1,  we  obtain  AE  =  12.7  kcal/mole  for  eq.  (4)  and 
AE  =  -5.2  kcal/mole  for  eq.  (5).  Since  we  find  6  to  be  slightly  more  stable  than  3  on  a  molecular 
level,  the  fact  that  3  is  the  form  that  is  observed  in  the  solid  phase  [4]  is  presumably  due  to 
stabilization  provided  by  crystal  lattice  interactions. 

In  seeking  to  understand  the  contrast  between  the  equilibria  in  eqs.  (4)  and  (5),  wherein  the 
former  favors  the  tetrazole  and  the  latter  the  azide,  it  is  useful  to  look  at  the  computed  electrostatic 
potentials  of  the  two  azides,  5  and  6;  these  are  shown,  on  the  molecular  surfaces,  in  Figures  4  and 
5.  The  minima  associated  with  the  nitrogens,  both  Vs,min(N)  and  VminCN),  are  listed  in  Table  3. 

An  interestins  feature  of  these  results  is  the  charge  distribution  in  the  azido  group.  Both 
Table  3  and  Figures  4  and  5  show  that  the  central  nitrogen  is  positive  in  character,  while  the  linking 
one  is  strongly  negative;  the  outermost  nitrogen  is  intermediate  between  these  extremes.  Thus,  in 
terms  of  the  resonance  description  of  the  azido  group  given  in  13A  -  13C,  the  dominant 
contributors  appear  to  be  13A  and  13B.  This  conclusion  is  consistent  with  the  structural  data  in 
Table  1,  as  well  as  with  earlier  interpretations  [32]. 

13  A  13B  13C 

Table  3  also  reveals  the  ring  nitrogens  in  5  to  be  more  negative  than  in  6.  In  this  respect, 
therefore,  our  finding  that  the  tetrazole  is  more  favored  energetically  in  eq.  (4)  than  in  eq.  (5)  is  in 
agreement  with  the  past  experience  mentioned  earlier  [5-9]. 

However  this  introduces  a  serious  question  concerning  4,  the  dinitro  derivative  of  3.  We 
did  verify  that  the  optimized  geometry  of  4  corresponds  to  a  trae  energy  minimum;  there  are  no 
imaginary  vibrational  frequencies  [33].  On  the  other  hand,  since  the  azido  tautomer  6  is  already 
shghtly  more  stable  than  the  tetrazole  3,  and  electron-withdrawing  groups  such  as  NO2  promote 
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the  formation  of  the  azide,  it  may  be  that 
energetically  over  4: 

0-,N  NO2 


the  azido  tautomer  14  is  considerably  favored 


4 


14 


(7) 


We  have  accordingly  computed  an  optimized  geometry  and  energy  for  14,  using  the  Gaussian 
92/DFT  procedure  described  earlier  [22-24].  The  resulting  energy,  -1000.37066  hartrees,  yields 
AE  =  -27.4  kcal/mole  for  eq.  (7).  This  confirms  that  the  introduction  of  nitro  groups  in  3  will  be 
expected  to  shift  the  tautomeric  equilibrium  significantly  toward  the  azide.  This  is  unfortunate 
from  the  standpoint  of  an  energetic  material,  because  of  the  well-known  proclivity  of  many  azides 
for  facile  decomposition  [32]. 

The  density  of  a  molecular  crystal  depends  upon  the  mass  and  size  of  the  molecule,  and 
upon  the  crystalline  structure  and  packing.  It  is  interesting  to  compare  1  -  3  to  phenanthrene  (15) 
and  acridine  (16).  The  molecular  weights  of  the  first  three  are  somewhat  less  than  those  of  the  last 
two,  yet  the  densities  of  1  -  3  are  more  than  30%  greater  than  those  of  15  and  16  [34].  Much  of 
this  difference  can  be  attributed  to  the  smaller  volumes  of  1  -  3  (Table  2),  all  of  which  are  at  least 
25%  less  than  the  201.4  K?>  of  15  and  199.3  of  16  [35]. 


This  is  fully  consistent  with  the  findings  of  Stine,  who  analyzed  the  crystal  structures  of 
more  than  2000  organic  compounds  to  obtain  effective  volumes  for  their  constituent  atoms  in 
various  bonding  environments  [36].  He  concluded  that,  in  general,  high  densities  are  promoted  by 
maximizing  the  numbers  of  nitrogens  and  oxygens  and  minimizing  the  numbers  of  carbons  and 
hydrogens. 

Molecular  interactions  are  also  expected  to  play  a  role  in  determining  the  densities  of 
organic  compounds.  We  have  found  to  be  an  effective  measure  of  this,  for  crystal  densities 
[25]  as  well  as  for  a  variety  of  other  macroscopic  properties,  including  critical  constants,  boiling 
points,  heats  of  vaporization,  partition  coefficients  and  supercritical  solubilities  [18-21,  25].  Since 
density  increases  with  [25],  the  large  magnimdes  of  the  latter  property  for  1  -  3,  greatly 
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exceeding  the  values  of  typical  organic  compounds  (as  pointed  out  above),  are  a  second  factor 
in  the  high  densities  of  1  -  3,  in  addition  to  their  relatively  small  sizes. 


Summary 

We  suggest  that  the  high  densities  found  crystallographically  for  1  -  3  can  be  attributed  to 
their  relatively  small  molecular  volumes  and  the  strong  intermolecular  attractions  arising  from  the 
highly  varying  electrostatic  potentials  on  their  molecular  surfaces.  3  and  4  can  be  involved  in 
tautomeric  equilibria  with  diazides,  and  in  the  gas  phase,  3  is  5.2  kcal/mole  more  stable  as  the 
diazide,  although  in  the  crystal  it  exists  only  as  the  tetrazole.  In  the  dinitro  derivative  of  3,  the 
energetic  preference  for  the  diazide  tautomer  is  greatly  enhanced,  indicating  that  its  usefulness  as  an 
energetic  material  is  open  to  question. 
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Figure  Captions 


Figure  1. 


Calculated  electrostatic  potential  on  the  molecular  surface  of  1.  Monochrome  r^ges,  in 
kcal/mple,  are;  dark  gray  for  V(r)  >  35;  light  gray  for  V(r)  from  35  to  0;  white  for 
negative  values  of  V(r). 


Figure  2. 


Calculated  electrostatic  potential  on 
kcal/mole,  are;  dark  gray  for  V(r)  > 
negative  values  of  V(r). 


the  molecular  surface  of  2.  Monochrome  ranges,  in 
35;  light  gray  for  V(r)  from  35  to  0;  white  for 


Figure  3. 


Calculated  electrostatic  potential  on  the  molecular  surface  of;  (a)  3  and  (b)  its_dmmo 
derivative  4.  Monochrome  ranges,  in  kcal/mole,  are;  dark  gray  for  V(r)  >  jo;  light 
gray  for  V(r)  from  35  to  0;  white  for  negative  values  of  V(r). 


Fiaure.  4.  Calculated  electrostauc  potendal  on  the  molecular  surface  of  a.  Monochrome  r^ges,  in 
°  kcal/mole,  are;  dark  gray  for  V(r)  >  15;  light  gray  for  V(r)  from  lo  to  0;  white  for 

negative  values  of  V(r). 


Figure.  5.  Calculated  electrostadc  potential  on  the  molecular  surface  of  6.  Monochrome  r^ges,  in 

"  kcal/mole,  are;  dark  gray  for  V(r)  >  15;  light  gray  for  V(r)  from  In  to  0;  white  for 

negative  values  of  V(r). 
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Table  1  DF/B-LYP/6-3 IG**  calculated  bond  distances  and  energies  and  experimentally-detennined 
_ bond  distances.^ _ _ _ _ _ 

calculated  bond  distances  (A)  experimentally-determined  bond  distances^  (A) 

_ and  energy  (hartrees) _ _ _ _ _ 


H 


H 


H 


X^.347 

N  N 

I  1 1.349 

C^423/C^ 

N  N  N 


£  =  -613.55104 


1.32^C>4.326 

N  N 

1.336|  |^•334 

1.322^C4.394  C^,333 

Nl397  N 

327]  ll.382  I  1.320 

Cu.32;^C^.334^G 

"1.3^N'^  N  1-32^j^ 


H 


1.318^C4.329 

N  N 

I.332j  l'-337 

1.34^C^  ^C^.343 

N  1.378'N  1.388  N 
.3141  ll-379  1  1.305 

C^1.32>C^.322^C 


H 


^.3^N' 


'N' 


H  H 

V/ 

.377  .N^ 


I  /b 

.N^./i.312 

^C'L473 

r\ 

H  H 


Nb 

\ 

N. 


H  H 
\  / 

1.340 


.N. 


N'  'C-33lV^^‘ 
I  1.328  1  Nb 

/l.3l5 

^C<460  Nc 

/\ 

H  H 


E  =  -592.60096 


H  H 

^C— 

/  \-395 

1.382 


U.391  / 

\^1_417^/ 


Nr 


336 


\ 


\  1.316 

-Nb 


Na  1-378 


H  H 

\  1.325  / 

c — c 

1.389/  \1.393 

l-250]vj^  N^-350 

^}C  \1.352 

1.299  V _ C  \  1.299 

Nw  /  \  ^Nb 

1.358  Na^'^*^ 


£  =  -591.39670 
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Table  1  DF/B-LYP/6-3 IG**  calculated  bond  distances  and  energies  and  experimentally-determined 
_ bond  distances  (continued).^ _ _ _ _ _ _ _ 

calculated  bond  distances  (A)  experimentally-determined  bond  distances^  (A) 

and  energy  (hartrees) _ _ _ _ _ _ _ 


0 

T 

0-N 

1.498'^ 


^1.241 

N-O 
1.369  /  1-238 


N 


Nb 


C— C 

/ 

N 


'N 


y.396  . 

\l_436/ 
/l.331 


N; 


•Na 

\  1.309 

^Nb 

1.380 


£  =  -1000.32701 


Nr 


H  H 

y  C  1-^21  IN  ^1.152 


1.285 


N, 


C  1-465 

/\ 

H  H 


N 


E  =  -592.58066 


\,.402  /H 

c— c 

/  \l-253 

N  N 

\  /l.342 

N— N—N 
c  b  a 


1.153 


N— N— N 

1.257 


£  =  -591.40491 


^The  experimental  data  for  1,  2  and  3  are  from  references  2  -  4,  respectively. 
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Table  2.  Calculated  and  experimental  properties  of  1  -  4.^ 


Molecule 

density 

surface 

area 

volume^’ 

— C _ 

^S,max 

Vs,min(N) 

^min(N) 

n 

9 

9 

^lot 

1 

1.689': 

163.9 

145.4 

52.2 

-43.2 

-78.1 

21.2 

181.7 

172.3 

354.0 

2 

1.715d 

162.0 

141.7 

53.7 

Na:  -38.9 
Nb;  -42.8 
Nc:  -28.5 

Na;  -72.0 
Nb:  -70.9 
Nc;  -52.3 

23.4 

111.9 

126.3 

238.2 

3 

1.674^ 

153.3 

133.7 

50.2 

Na;  -25.4 
Nb:  -41.4 
Nc:  -44.0 

Na:  -47.7 
Nb;  -68.8 
Nc:  -73.3 

24.7 

130.6 

137.0 

267.6 

4 

— 

200.0 

174.3 

72.0 

NV  -18.6 
NV  -31.3 
Nc:  -32.2 

Na:-48.6 
Nb;-54.9 
Nc:  -57.9 

21.6 

358.6 

o 

64.3 

422.9 

... _ _ _ — - - - - - -  9  2  2 

^Units',  density,  g/cm^;  area,  A^;  volume,  A^;  all  potentials  and  FI,  kcal/mole,  <J_  and 


(kcal/mole)2.  jhe  labeling  of  the  nitrogens  follows  Table  1. 
l^Reference  36. 

'^Reference  2. 

^Reference  3. 

^Reference  4. 
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Table  3.  Most  negative  electrostatic  potentials,  Vs,min(N)  and  VminlN),  associated  with  the 


nitrogens  in  the  azides  5  and  6.^ 

Molecule 

Vs.min(N),  kcal/mole 

Vmin(N),  kcal/mole 

C 

Nring* 

-32.6 

Nring- 

-122 

Na: 

-29.3 

Na 

-57.9 

Nb: 

none  found 

Nb 

none  found 

Nc: 

-19.7 

Nc 

-28.0 

6 

Nring: 

Na: 

-28.6 

-49.2^ 

Nring 

Na 

-65.9 

:  -68.3 

Nb: 

none  found 

Nb 

:  none  found 

Nc: 

-19.9 

Nc 

;  -28.1 

aThe  labeling  of  the  nitrogens  follows  Table  I . 

bThe  negative  surface  region  associated  with  each  Na  overlaps  that  of  the  other,  and  there  is  only  one  loc 


minimum. 
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Fi<^ure  3  Calculated  electrostatic  potential  on  the  molecular  surface  of:  (a)  3  and  (b)  its  dinitro 
®  derivative  4.  Monochrome  ranges,  in  kcal/mole,  are:  dark  gray  for  V (r)  >  35;  tight 

gray  for  V(r)  from  35  to  0;  white  for  negative  values  of  V(r). 


Figure.  4. 


Figure. 


Calculated  electrostatic  potential  on  the  molecular  surface  of  5.  Monochrome 
ranges,  in  kcaFmole,  are:  dark  gray  for  V(r)  >  15;  light  gray  for  V(r)  from  15  to 
0;  white  for  negative  values  of  V (r). 


Calculated  electrostatic  potential  on  the  molecular  surface  of  6.  Monochrome 
'  ranges,  in  kcal/mole,  are:  dark  gray  for  V(r)  >  15;  tight  gray  for  V(r)  from  15  to 
0;  white  for  negative  values  of  V(r). 


